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We study a planar four-Fermi Gross-Neveu model in the presence of a tilted magnetic field, with 
components parallel and perpendicular to the system's plane. We determine how this combination of 
magnetic field components, when applied simultaneously, affects the phase diagram of the model. It 
is shown that each component of the magnetic field causes a competing effect on the chiral symmetry 
in these fermionic systems. While the perpendicular component of the magnetic field tends to make 
the chiral symmetry breaking to become stronger, the effect of the parallel component of the field in 
these planar systems is to weaken the chiral symmetry. We show that this competing effect, when 
combined also with temperature and chemical potential, can lead to a rich phase diagram, with 
the emergence of multiple critical points and reentrant phase transitions. We also study how the 
presence of these multiple critical points and reentrant phases can manifest in the quantum Hall 
effect. Our results provide a possible way to probe experimentally chiral symmetry breaking and the 
corresponding emergence of a gap (i.e., the presence of a nonvanishing chiral vacuum expectation 
value) in planar condensed matter systems of current interest. 

PACS numbers: 11.10.Kk,71.30.+h,11.30.Qc 



I. INTRODUCTION 

The interaction of fcrmions with an external magnetic field is expected to be associated with phenomena like 
metal to insulating phase transitions in semiconductors quantum Hall effect and transport properties in 
superconductors [3], just to mention a few examples in the context of condensed matter physics. Among these 
systems, planar ones are of particular interest and they include, for example, high-temperature superconductors, 
organic thin films and most recently, the physics of graphene [|| . Most of these systems have an excitation spectrum 
that can be well described by relativistic Dirac-like fermions. A typical example is the physics of graphene, in which 
electron transport is essentially governed by the Dirac's relativistic equation. It is then expected that these type of 
planar systems can be appropriately described by quantum field theory models in two spatial dimensions. In this 
context, Gross-Neveu (GN) type of models Q in 2 + 1-dimcnsions are very popular not only because of their simplicity, 
but also for their ability to capture many of the relevant physics exhibited by planar fermionic systems in general. 
These models can have both a discrete, ip — > Jsip, or a continuous, ip — > eKjp(ia"f^)ijj , chiral symmetry. They consist 
of N massless four-component fermions self-interacting through a local four-Fermi interaction. 

One of the most basic motivation for studying the properties of fermionic planar systems is whether we can have a 
metal to insulating type of transition under the variation of external parameters, like temperature, chemical potential 
and external fields, e.g., an external magnetic field. The metal-insulator transition is directly connected to the 
presence of a non-vanishing gap in the excitation spectrum for the fermions and how this gap eventually vanishes 
in the presence of external conditions. In the quantum field theory context, a mass term for the fermions can be 
generated dynamically in the absence of external conditions, like in the original GN model, and it is associated with 
the breaking of the chiral symmetry in the model. 
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FIG. 1: A schematic representation of the system under study, consisting of fermions in a plane and in the presence of a 
magnetic field B that makes an angle (f> with respect to the plane. 



In this work we make use of the original GN model with discrete chiral symmetry and study the chiral phase 
transition when the system is in the presence of an external tilted magnetic field, i.e., a magnetic field that includes 
both parallel and perpendicular components that are being applied simultaneously to the system (see Fig. [T]). This 
model has been previously extensively analyzed when in the presence of a magnetic field applied perpendicular to 
the plane (see, for example, Refs. [f| 0] for ones of the earlier studies) and, more recently, in the presence of an 
in-plane magnetic field P. But, surprisingly, a study of both type of fields, parallel and perpendicular to the system's 
plane and applied simultaneously is still lacking. Perhaps, the reason for this negligence in the literature was because 
not much was expected for the effects on the system symmetry due to a parallel component of the field. However, 
there are good reasons for making this study. It is well known that a perpendicular magnetic field B±_ leads to an 
enhancement of the chiral symmetry broken region (thus, it tends to increase a gap in the excitation spectrum). 
This is the magnetic catalysis phenomenon that can occur for all values of the four- Fermi interaction strength 0, 0] . 
However, an in-plane magnetic field B\\ has been shown to produce an opposite effect [H, leading to a reduction of the 
chiral symmetry broken region. In other words, a parallel component of the magnetic field tends to act on restoring 
the chiral symmetry for large enough values of B\\ . This is due to the fact that for a planar fermionic system, the 
in-plane magnetic field contributes like a Zeeman energy term, leading to an effective polarization of the system. The 
polarized planar fermionic system tends to exhibit less chiral symmetry. 

Since the application of perpendicular and parallel (in-plane) magnetic fields leads to opposite effects on the chiral 
symmetry, it is reasonable to expect that the combination of both components, when applied simultaneously to the 
system, to possibly lead to new patterns for the chiral phase transition. We will show in this work that this indeed 
can happen under appropriate choice of parameters. In this paper we show that the application of both B± and £?ii, 
with the combined effects of temperature T and chemical potential (j,, can lead to a phase diagram that can display 
a rich structure, with the possibility of exhibiting multiple critical points and reentrant phases. We show that this is 
a consequence of the competing effects on the system's chiral symmetry coming from each component of the applied 
oblique external magnetic field. 

This paper is organized as follows. In Sec. |H] we briefly revise the four-Fermi GN model used in this work. In 
Sec. Mil the thermodynamic effective potential for the fermionic scalar condensate (ipip) is derived in the mean- field 
theory approximation and results for the phase diagram for the model are displayed. In Sec. HVI we study in particular 
the quantum phase transitions in the model and derive the Hall conductivity. We also discuss the impact that the 
presence of the possible multiple critical points can have on the Hall conductivity, pointing out to a possible way for 
characterizing experimentally the presence of quantum critical points. Finally, in Sec. |Vj we give our conclusions and 
discuss possible applications and extensions of this work. 



II. THE FOUR-FERMI MODEL IN THE PRESENCE OF AN EXTERNAL MAGNETIC FIELD 

We work with a four-Fermi model, in 2 + f-dimcnsions, of the GN type and with N flavors, described by the 
Lagrangian density [|| 



£ =EE 



2N 



(2.1) 



where A is the coupling and $ = r f J d v , v = 0, 1,2, with the gamma matrices being 4x4 matrices, where here we 
follow the representation given, e.g., in Ref. @ for fermions in 2+1-dimensions. The model (|2.1[) possesses a discrete 
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chiral symmetry, tp — > 75^, — > —'075, with the 75 matrix defined as in Ref. Q. This discrete chiral symmetry is the 
one considered in this paper, along with its breaking and restoration conditions. The chiral symmetry is broken when 
a gap (a mass term for the fermions) is generated. Since the gap corresponds to a nonvanishing vacuum expectation 
value for the chiral operator, {ifrip), it is then convenient to rewrite Eq. (|2.ip in terms of an auxiliary scalar field a, 



L[^,a}= £ ^(^-<r)^--<r 2 , (2.2) 

s =t,J.i=l 

and where a and the chiral operator are related, from the saddle point solution for a, by a = —(\/N)ij)ip, with implicit 
sums over the spin and flavor. In terms of the Lagrangian density (|2.2[) . the discrete chiral symmetry is now expressed 
as ip — > 7505 "0 — * -075 a- n d o — > — a. 

With the Lagrangian density Eq. (|2.2j) expressed in terms of the auxiliary scalar field, it becomes relatively simpler to 
study the phase diagram of the model, which can be made through the study of the behavior of the vacuum expectation 
value of a as a function of the temperature, chemical potential and the external magnetic field. This study can be 
made by evaluating the effective potential, or Landau's free energy density, for a constant (background) configuration 
for the auxiliary scalar field, V e e(a c ). The effective potential can be obtained by integrating out the fermion fields 
and fluctuations around the scalar background field Q. From the effective potential, all the thermodynamics for the 
model can be derived. In particular, a chiral phase transition is signaled by a non- vanishing vacuum expectation value 
for a, ij c = (a), which is a (global) minimum of V c s{a c ). Both temperature and chemical potential are introduced 
through the grand canonical partition function definition, 



Z {P, =Tre X p I- j3 (H-fiQ)} , (2.3) 

where /3 is the inverse of the temperature 1 , /j, is the chemical potential, H is the Hamiltonian corresponding to Eq. (|2.2[) 
and Q = J d 2 xtpjoip is the conserved charge. Furthermore, we can also study the effect of an external magnetic field 
applied to the system by coupling it to the fermions in (|2.2[) through the usual way. For an arbitrary external magnetic 
field that can have components both parallel and perpendicular to the system's plane, each component will couple to 
the fermions differently. The perpendicular component of the field couples to the fermions orbital motion, generating 
the Landau levels for fermions in a magnetic field @, 0] ■ The parallel (in-plane) component of the field couples only 
to the spins of the fermions, producing an intrinsic Zeeman effect fl(]| . The effective potential, or free energy density, 
is defined as usual from the grand canonical partition function (|2 .3[) by 



V eS = ~\aZ, (2.4) 

where V is the volume (or, more specifically, the area in the present problem of a planar system). 

Transforming Eq. (|2.3[) in the form of a path integral, in the imaginary-time (Euclidean) formalism of finite tem- 
perature field theory , one obtains 



Z= DaY[D^Dip exp{-S E [ip,il>,a]} , (2.5) 

J s 

where the functional integration is performed over the fermion fields satisfying the anti-periodic boundary condition 
in Euclidean time: rpj (x, r) = —tpj (x, r + /?) , while over the a it satisfies periodic boundary conditions. The Euclidean 
action Se[4>, ip] m Eq. (|2.5j) . from the Lagrangian density Eq. (|2.2j) in the presence of an external magnetic field with 
components B± and Bu , perpendicular and parallel to the plane, respectively, is given by 2 



1 Along this work, we assume the natural units, kg = h = c = 1. 

2 Note that in condensed matter systems there should also be a Fermi velocity vp factor multiplying the space components of the gamma 
matrices and the coupling constant A in this expression Q (for example, in graphene vp rj c/300). vp is omitted in all our expressions, 
but we will re-established it when making estimates in the context of condensed matter systems later on. 
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where A x and A y are the external electromagnetic vector potential components and a s jo g [1bB\\/2 is the corresponding 
Zeeman energy term, with cr-|- = 1, <rj, = — 1, g is the spectroscopic Landc factor and \ib is the Bohr magneton. By 
choosing a gauge where the three-dimensional vector potential is given, for example, by A — (0, B±x, B^y), we see 
from Eq. (|2.6p that B± couples only to the orbital motion of the fcrmions and it will result in the Landau levels for 
the system in this magnetic field. The parallel magnetic field couples to the electrons' spin and produces the Zeeman 
energy term in Eq. (|2.6[) . From the form of the Zeeman energy term in Eq. Q2.6p . we see that it can be added to the 
chemical potential. Thus, in the action Se, we can define an effective chemical potential term of the form, 



N N 

E E rt ? Tov,- = E E 0* + y 9»bb^ 

= + V , (2.7) 

where = fx + Sp,, and /Ltj, = fx — S/j,, with 8pL = g pbB\\/2. The effect of the in-plane component of the magnetic field 
is then felt through the asymmetry term 5fi in the chemical potential for spin up and down, resulting in a polarization 
of the system. This polarization produced by the Zeeman energy term is directly proportional to Bu . Note that in the 
language of condensed matter systems, the in-plane magnetic field produces the analogous of an asymmetrical doping, 
or imbalance between the chemical potentials of the fermions (electrons) with the two possible spin orientations. This 
asymmetry can be produced by a doping process that alters the densities of spin up and spin down. Physically, the 
chemical potential [i can be interpreted as to account for the extra density of electrons that is supplied to the system 
by the dopants, while Sfi (produced by the in-plane magnetic field component) measures the amount of asymmetry 
introduced. 



III. EFFECTIVE POTENTIAL AND PHASE STRUCTURE 



All the thermodynamics, magnetic properties and phase structure of the model (|2.2I) . as a function of temperature, 
chemical potential and magnetic field, can be determined from the effective potential (|2.4[) . It is useful to start by 
briefly reviewing the derivation of the effective potential and its main results for the GN model, when including first 
only temperature and chemical potential, by then, later, by also including the external magnetic field. 



A. The effective potential and phase structure at finite T and n 



From the grand canonical partition function Eq. (|2.5[) , with the action (|2.6[) (taking initially in the absence of 
an external magnetic field), for a constant background scalar field (a) = a c , in the mean-field approximation, or 
equivalently, from the leading term in a 1/N expansion, or the large- N approximation [l~2l Il3j . we obtain that the 
effective potential V c s for er c , at finite T and n, is given by 



N 



V eS (a c ,T,fi) = — a 2 c -2NT £ 



d 2 p 
(2tt)2 



In 



(w n - + p 2 + <y\ 



(3.1) 



where uj n = (2n+l)7rT are the Matsubara frequencies for fermions. Performing the sum over the Matsubara frequencies 
in Eq. ([5TT]) . we find 



y off (a c ,T, M ) = ^ra 2 c -2N 



2A 

2iVT 



d 2 p 
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(27T) 2 



(27T) 
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2 E P 



(3.2) 
(3.3) 
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where E ± = E p ± \i and E p = \J\> 2 + a 2 . 

At zero temperature and chemical potential the effective potential becomes 



V cS {o c ) = ^a 2 c -2N 



(3.4) 



The momentum integral in Eq. (|3 .4[) is ultraviolet divergent. By using a momentum cutoff A to regulate the 
divergent term in V e s(a c ), we can define a renormalization condition for the coupling as 



1 d 2 V eS {a c ) 



X R (m) N da 2 

where to is a regularization scale. Equation (|3.5[) gives (neglecting terms of 0(1/A)), 



(3.5) 



1 1 A 2m 

Ar(to) A 7T 7T 



(3.6) 



Next, by defining a renormalization invariant coupling as 14 j 



1 



1 



2to 



Xr A_r(to) 7T 



(3.7) 



we obtain that in terms of A/j the renormalized effective potential reads (after subtracting an irrelevant field indepen- 
dent divergent vacuum term) 



V eS , R (a c ) = —a 2 + — \a c \ 3 



(3.8) 



The minimum a c of ~V e fj,R{pc) is given by 



e>o , X R < 0, 
, Xr > , 



(3.9) 



where <7o = 7r/|A#|. Thus, dynamical chiral symmetry breaking, in the absence of external fields, only occurs for a 
negative coupling constant. 

Using (|3.8p and performing the (finite) momentum integrals for the temperature and chemical potential dependent 
terms in Eq. (|3.3j) . the renormalized effective potential at finite T and /i becomes 





Ksff,R.(Cc) 
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7T I 



,-/8(<7 c -a.) 



]+Li 3 [-e 



-/3(<x c +a0 
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where Li„(z) is the polylogarithm function and it is defined (for v > 0) as [1£ 



(3.10) 



fe=i 



From Eq. (|3.10|) . we can derive the gap equation, which is defined from 

d 



da. 



V e s(cr c ,T, (J,) 



a c =a c (T,ti) 



= o, 



(3.11) 
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which gives 



G C = <7 C (T = 0,(1 = 0) 



In 



1 



111 



1 



-/9(ffc-M) 



(3.12) 



The behavior of a c is well known 

[UEl. At T = 0, there is a critical chemical potential fx c at which the chiral 
symmetry is restored (we are considering the case here where chiral symmetry breaking happens, i.e., for A/? < 0). 
This can be shown to happen through a first order phase transition, defined by the relation T4ff i R,(cr c = 0, fi c ) = 
Kff,R(c c = ero,/j c ). This yields the critical chemical potential /x c = uq. Thus, at zero temperature, we have the 
ground state of the system given by 



fcxo, for M < Mc , 

[ 0, for /z > fi c . x ' 

At finite temperature and chemical potential, there is a critical curve (T C (T, (i) = 0, which is obtained from Eq. (|3.12l) . 
This gives a line of second order phase transition in the (//, T)-plane, starting at the critical point (ii = 0, T = T c ), 
where T c = oo/(21n2), and ending in a first order critical point at (u = /i c ,T = 0). Note that when going beyond 
the mean-field approximation this structure can slightly change [17] . leading to a first-order critical line merging 
with the second-order critical line at a tricritical point. In this work, we will restrict our analysis at the mean-field 
approximation only. 



B. The effective potential for the GN model at finite T, fi and in a tilted magnetic field 



Let us now investigate the effective potential for the system in the presence of a tilted magnetic field, with com- 
ponents (Bn,Bj_). The effective potential can once again be derived from Eq. (|2.5p and it can also follow directly 
from Eq. (|3.ip or (|3.3p . In the case the magnetic field has only an in-plane component, B\\ ^ and B± = 0, we have 
a shift of the chemical potential and we need to distinguish in the effective potential, besides the contribution from 
particles and anti-particles, also the contributions from particles with spin up and spin down, with chemical potentials 
m = (i + 5(i and /Uj. = jx — 5jjL, respectively, with <5/i = g hbB\\/2. Along this work we will be considering (i and 5(i as 
positive quantities. Equation (|3.3[) for the renormalized effective potential, in this case, becomes 



V eff (o- c ,T, (j,f, m) 



N 2 N ! ,3 

2^ + 3^l 
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(3.14) 



where E. 



T4 



E p ± |^,^|. A detailed analysis of the phase structure in this case was performed in Ref. 



The 



main effect of the in-plane magnetic field is to further contribute towards chiral symmetry restoration. In this case, 
it can be shown that the chiral symmetry will remain always restored (at any T and (i) for a critical asymmetry 
S(i c > (i c = do- 
Finally, we can also include in the above equation (|3.14[) a perpendicular component for the magnetic field, B±^ 0, 
by properly accounting for the Landau energy levels and degeneracies. For instance, in Eq. (|3.14p . we can take @!0 
E p = \/p 2 + &c \f^k\eB±_\ + a 2 , k = 0,1,2,..., the momentum integrals are replaced by a sum over the Landau 
levels, with a density of states \eB±\/ (2ir) and by also accounting for the degeneracy of the k > 1 Landau levels, we 
obtain 
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V eSiR (a c ,T,n,B ± ,5n) = — C oi u i + 1 ~ 



TV cn,nA c ' x ' ^ 2X R tt S V 2'2|efl±| 7 2ir 

\eB 



] {in (l + e~^-^) + In (l + e -«l"<>l+/t>) 
2 £ In (l + (V^\eB ± \- H ) 



k=l 



where C(s,a) is the Hurwitz zeta function [Ta |. 



2]f> (l + e-KV^+ 2fe l^l+^)^ 
fc=i ^ ' 
(Mt->N)}. (3-15) 

oo 1 



=0 



Equation (|3.15 1 can be seen as a generalization of the mean- field theory (or large- N approximation) result derived, 



e.g., in Refs fA |14||. for the case considered here, where we have a parallel component of the magnetic field (or 
equivalently, an asymmetry between the spin up and spin down of the fermion field). 

An important feature concerning the 2+1-dimensional GN model in a magnetic field, when B± ^ and taking 
B\\ = 0, is that it can be shown that a nonvanishing global minimum can always develop for the effective potential, 
even when the couplingin (|3.15|) is positive. This is the phenomenon of magnetic catalysis, found and studied at 
length, e.g., in Refs. [si?}. The behavior for the global minimum a c for small values of B±. at Bn = T = /i = 0, has 
been determined to be given by @ 

fs?|Afi||&B±| + ." , A fi >0, 



i i \eB 



12a; 



<3 ' 17) 



Thus, the component of the magnetic field perpendicular to the plane always tend to break (or enhance) the chiral 
symmetry. On the other hand, as we have already commented above, a parallel (in-plane) magnetic field component 
Bn 7^ has been shown Q to always tend to restore the chiral symmetry. Therefore, when both components of the 
magnetic field are present, one can expect that there must be a competing effect between the in-plane component 
Bu and the perpendicular component B±. While the former tends to suppress the chiral symmetry, the latter 
tends to enhance it. We then expect that the interplay of the two components of the magnetic field, when applied 
simultaneously and independently to the system, will possibly lead to a rich structure for the phase diagram. This 
indeed will be true, as we will show below. As we will show, the simultaneous and independent application of external 
magnetic fields in parallel and perpendicular to the system's plane can generate multicritical points and reentrant 
phases. These are results not seen in those cases when only one component of the magnetic field is present. It is 
quite surprising that such study of the simultaneous effect of both components of the magnetic field in the GN in 
2+1-dimcnsions has not been done before, despite the fact that this is a kind of physical situation that has been of 
increasingly interest in the context of experiments with planar type of condensed matter systems in the laboratory. 
We will comment more on this in the conclusions section. 

C. The phase diagram in the presence of a tilted magnetic field 

We now turn to the study of the phase diagram for the 2+1-dimensional GN model, with the effective potential 
given by Eq. (|3.15l) . The equation for the gap a c , which is determined from the equation defining the minimum of 
the effective potential, 



-J— V e S,B.((Tc, T, [1, Bj_, 5(1) 



, (3.18) 
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and that generalizes Eq. (13 . 12[) . becomes 



V2 
2ir 



+ 



47T 



1 



1 



2' 2|eB_L| 
1 



2tt 



-mt) + 1 e Pi<Tc+IH) _|_ i 



+2cr c ]T 



- 1 (^ + 2fc|eS ± | 



/(V*c+2fe|eBx|+Mt) 



lJ 



e ^(\/*2+2fe|eSx|-Mt) 



(3.19) 



A second-order phase transition critical point is found when, by changing any of the external parameters T, /i, Bj_ 
or 5/i(_B||), (j c changes continuously and vanishes at the critical point. For a first-order phase transition, we have, 
instead, that a c is discontinuous at the critical point. We are interested in obtaining the critical points by varying 
a set of parameters, while others are kept fixed. For example, we can initially fix the values for the in-plane and 
perpendicular to the plane components of the magnetic field, while varying the temperature and chemical potential. 
For fixed values of B± and S[x (i.e., for a second-order transition line in the (/x, T)-plane then follows from those 
critical points for which a c vanishes in a continuous way along this critical line. In a first-order transition, the effective 
potential develops different minima, ^ <7c , where one of them is a local minimum, while the other is a global 
minimum. These minima can get degenerate for some values of the parameters. The first-order transition line is then 
determined by the condition of degeneracy of the minima of the effective potential, 



V c f f , R (ai 1 \T c , f i c ,B ± ,S^) = V eS , n {af\T c ,ii c ,B x ,5^) . (3.20) 

In general, one of the minima is the trivial solution, a c — 0, though this may not be always true, depending on the 
parameters chosen in Eq. (|3.15p . as we will verify explicitly below. In all cases, the critical lines are determined 
numerically, by solving the gap equation (|3.19p and by also verifying the condition p.20j) . The point where the 
second-order critical line meets the first-order critical line defines a tricritical point in the phase diagram. 




FIG. 2: The phase diagram in the {fi, T)-plane for Xr < 0. The plot on the left shows the critical curves in the case of B± — 
and for different values of <5/x. The plot on the right shows the critical curves for the case of By = (Sfj, = 0) and for different 
values of B±. The solid lines are second-order critical lines and the dotted lines are first-order critical lines. The tricritical 
points are marked by the large dots in the plot on the right. All quantities are expressed in units of ao — 7t/|Ak|. 

Let us initially show the effect on the phase structure of the system when only one of the components of the 
magnetic field is applied to the system, i.e., for the case of B» ^ 0, B± = (corresponding to a tilted angle of zero 
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degrees of B with respect to the system's plane) and for Bn — 0, B± ^ (corresponding to a tilted angle of ninety 

degrees for B). In Fig. [5] we show the effect of each component of the external magnetic field, when applied to the 
system separately, in terms of the phase diagram in the (fi, T)-plane. In the two plots shown in the figure, the external 
region to the curves is a chiral symmetry restored phase region, a c = 0, while the internal region to the curves is a 
chiral symmetry broken phase region, a c 7^ 0. For simplicity, we are not showing the mctastablc regions that would 
appear around the first-order transition lines in the plot on the right. These metastable regions will be discussed and 
shown below. Figure [Ha) is for the case where the perpendicular component of the magnetic field is absent, B± = 0. 
In this case, the critical lines for different values for the asymmetry 5[i (Bu 7^ 0) are found to be all of second-order. 
The only first-order transition point happens at T = and (i c = <jq. For a critical value of the asymmetry <5/i c = cro, 
we find that a c = for any value of the temperature. Note that this critical value for the asymmetry has the same 
value as for fi c , found in the absence of the Zeeman energy term fl6l ] . This is a consequence of the symmetry between 
the Zeeman energy term (Sfi) and the chemical potential. Figure EJb) is for the case where the in-plane component 
of the magnetic field is absent, B\\ = 0, but for B± 7^ 0. The structure obtained is the one previously found in the 
earlier papers 0, 0], where, besides the second order transition line, now a first order transition line also appears, 
with a tricritical point between them. The effect due to magnetic catalysis is clear, showing that the larger is B±, 
the larger is the symmetry broken region. It is then clear from the results shown in Fig. [2] that the in-plane magnetic 
field tends to act to restore the chiral symmetry, while the perpendicular magnetic field tends to enlarge the chiral 
symmetry broken region. 
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FIG. 3: The phase diagram in the (/x, T)-plane for Ar < 0. The plot on the left shows the critical curves in the case of fixed 
je_Bx| = 7<To and for different values of Sfi. The plot on the right shows the critical curves for the case of fixed Sfi — 1.5<ro and 
for different values of |ei?x|. The solid lines are second-order critical lines and the dotted lines are first-order critical lines. The 
tricritical points are marked by the large dots in the plots. 

In Fig. [3] we show the effects of increasing one of the components of the magnetic field, while the other component 
is kept fixed. In Fig. EJa) the perpendicular magnetic field is kept constant at the value of |eBj_| = 7<7q and the 
in-plane magnetic field (or cquivalcntly, the asymmetry Sfi) is varied. In Fig. EJb) we show the opposite case, where 
we kept fixed the in-plane magnetic field, with resulting asymmetry given by Sfi — 1.5ao, while the perpendicular 
component of the magnetic field is varied. A noticeable feature that we now observe in each of these cases, not seen 
in those cases where only one of the components of the magnetic field is present, is the possibility of emergence of 
multiple critical points and reentrant phases. In the case of fixed B±, Fig. [3Ja), after some critical value for the 
asymmetry, which for the parameters taken for this plot is 8(1 ~ 0.895cto, a second tricritical point emerges with 
an associated first-order transition line. In the plot shown in Fig. I^b), we have purposivcly considered a value of 
asymmetry of 5/i = 1.5<7q- Recall that for an asymmetry above the critical value of 5/i c = o$ and when B± = 0, the 
chiral symmetry would always be restored However, as we start increasing B±, a symmetry broken region starts 
to emerge after some large enough value for the perpendicular magnetic field. This region of chiral symmetry broken 
phase that appears has also the form of the multicritical curves shown in Fig. [3ja). For a sufficiently large value of 
B±, which for the parameters used in the plot of Fig. [3Jb) is \eB±\ ~ 35ctq, one of the additional tricritical points 
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FIG. 4: The phase diagram (including the metastable lines) in the (/i, T)-plane. The perpendicular magnetic field is fixed in 
|eBj_| = 7o"o- The panel (a) on the top left shows the critical curves in the case of Sfi = 0. The panel (b) on the top right 
shows the critical curves for the case of Sfj, = 0.8<ro- The panel (c) on the bottom left shows the critical curves for the case 
of Sfi = 0.897oo. The panel (d) on the bottom right shows the critical curves for the case of 8fi — 1.5cto- The solid lines are 
second-order critical lines, the dotted lines are first-order critical lines and the dashed lines are the metastable critical lines. 
The tricritical points are marked by the large dots in the plots. It is also shown the schematic shape of the effective potential 
in each region of the phase diagram. 



(and associated first-order transition line) disappears and the phase diagram shows again a standard form with a 
second-order critical line, a first-order transition line and a tricritical point in between them, a structure just like the 
ones shown in Fig.[2Jb). The combined effect of both components of the magnetic field, parallel and perpendicular to 
the plane, is able to lead to reentrant phase transitions, as seen from the results in Fig. [3J whenever more than one 
tricritical point is present. This feature will be even more apparent when we consider the zero temperature limit and 
look at the possibility of quantum phase transitions at fixed values of B± and B\\ and varying chemical potential, a 
situation we will analyze in more details below, in Sec. H*Vl 

For completeness, in Fig. U we present the detailed structure for some of the phase diagrams shown in Fig. [3J 
including now the metastable lines. The metastable lines give the points for which a local minimum of the effective 
potential first appears or disappears. These lines are the ones shown on the left and on the right of the first-order 
critical line. Again, from the results shown in Figs.^Jc) andUJd), we can see the presence of reentrant phase transitions 
at some fixed temperature and varying chemical potential (e.g., the analogous of varying the doping in condensed 
matter systems). Reentrant phase transitions are always associated with more than one critical point. Multiple 
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critical points can appear with more frequency as we decrease the temperature and, in particular, as we enter the 
quantum regime (T = 0). In the next section we will analyze this situation more carefully. 

IV. QUANTUM PHASE TRANSITIONS PATTERNS 

In order to better understand the reentrant phases found in the previous section, it is convenient to first look at 
those cases where T = 0, i.e., for the quantum regime. Taking the T = limit in Eq. p,15[) . we obtain 



-^V c ff,R(a c ,T,fi,B±,5fi) T ^ 



al V2\eB ± \W ( 1 o\ 
2X R 7T ^ I 2'2\eB A 



2tt 



(A*t - K\) o (Mt - kcl) + (Kl - W\) e - kd) 



47T 

fc=i 

+ 2 >T - y/2k\eB ± \ + ai) 6 (ft - a 2 c ) 

k=l 



(4.1) 



where 



^'jnav Int 



2|eS. I 



(4.2) 



and Int(x) means the integer part of x. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
yU/CT 

FIG. 5: The phase diagram in the (fx, <5/i)-plane and at T = 0. The perpendicular magnetic field component was fixed at the 
value of |e£>x| = <Jq. The different shadows of gray indicate decreasing values for the vacuum expectation value a c . Going 
from dark to light shadows, we have a c /ao — 1.0617, 0.8144, 0.3669, 0.2883, respectively. Outside the shadow areas a c = 0. 
All solid lines indicate first-order phase transition lines. 



We can now look at the phase diagram in the (fj,, (S^)-plane at fixed B±. An example is shown in Fig. [5j Note from 
Fig. [5] that for fixed values for the asymmetry S[i (i.e., -By), we can find multiple reentrant phases. Examples of how 
these reentrant phases manifest from the effective potential are shown in Fig. [S] for two specific values of asymmetry, 
5fi = 0.7<7o and 5/i = 0.9cto, respectively. 
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FIG. 6: The effective potential for a fixed value of |eBj_| = cr and for different values of chemical potential fj, and asymmetry 
Sfi. The plot on the left is for 5fx = 0.7cro and for fi/ao — 0.2467 (solid curve), fi/ao ~ 0.8633 (dashed curve) and fi/ao ~ 0.9845 
(dotted curve). For fi/ao > 0.9845, there is only one global minimum at a c = and chiral symmetry is restored. The plot on 
the right is for Sfi — 0.9<to and for jti/uo = (solid curve), (j,/ao — 0.3172 (dashed curve), fi/ao — 0.6633 (dash-dotted curve) 
and fi/ &o — 1-1282 (dotted curve). For /i/ao < 0.3172 and for fx/ao > 1.1282, there is only one global minimum at a c = and 
chiral symmetry is restored. 

From the forms of the effective potential shown in Fig. [51 we see that two forms of (first-order) reentrant phase 
transition are possible at some value for the asymmetry and when the chemical potential is increased. We can have 
a reentrant phase transition starting at a chiral broken state, a c = <Ji, go to one or more intermediate chiral broken 
states with a c = (72, &i > 02, and then to a chiral symmetric state, a c = 0. This is the case shown, e.g., in the first 
panel in Fig. [6] The other possibility is, starting from a chiral symmetric state, a c = 0, go to one or more intermediate 
chiral broken states, with a c ^ 0, and then back to the chiral symmetric state. This is, e.g., represented in the second 
panel in Fig. [6] 

A. Quantum phase transitions in the fully polarized system 
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FIG. 7: The values assumed for the vacuum expectation value a c along the line 8fi — fi in Fig. [5] 

It is interestingly to notice from Fig.[S]that the reentrant phases are maximal (occur more frequently) when /i = <5/i. 
This corresponds to the full polarization regime, ^ = 2fi and fj,^ = 0. In Fig. [7] we show the values assumed for the 
vacuum expectation value a c along this full polarization regime. By increasing the chemical potential, the value of a c 
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can decrease discontinuously through many intermediate first-order phase transitions. This type of pattern is found 
to be exclusive of the full polarization regime. 



1.2 r 



1.0 




FIG. 8: The values assumed for the vacuum expectation value a c along the line <5/x = fj,, when |eBx| = a% and T = 0.06cro. 

Though this pattern of multiple phase transitions can happen at any nonzero value of B± in principle, this structure 
is not stable when thermal fluctuations are accounted for. Leaving the quantum regime T = and by increasing 
the temperature, the intermediate transitions going from a value of a c = o\ to a value (T2 (where a\ > 02) quickly 
disappears. This can be easily understood from the results shown for the effective potential in Fig. [51 These interme- 
diate first-order transitions are in general characterized by a small potential barrier. By increasing the temperature, 
these intermediate minima in the potential are quickly smooth out. The smaller is Bj_, the smaller is the temperature 
needed to smooth out the intermediate transitions. For example, in Fig. 13 we have the case of \eB±_\ = <r 2 and for a 
temperature of T = O.O60O7 which are taken as representative values. All intermediate transitions are smoothed out 
and a c becomes mostly a continuously varying function, remaining only the strongest first-order transition, which is 
the first one shown in Fig. [7] and that also remains in Fig. [8l By further increasing the temperature, the transition 
becomes just a second-order one, with a c continuously varying till vanishing completely at a given critical value of 
chemical potential. This corresponds to crossing one of the lines of second-order phase transition. 

From the expression for the effective potential in the T = limit, Eq. (|4.1j) . we can find approximate analytical 
expressions for the vacuum expectation value d c in the large magnetic field regime, with |e-Bj_| 3> \J? and |eBj_| 3> <J 2 a - 
For instance, in the /x = Sfz = case, we find, valid for either positive or negative values of Xr, 

M r = o,^ = o,, = o),^_ 7 ^L_, (4 „ 

where £(1/2) ~ —1.46035, while for the case when there is only the perpendicular component of the magnetic field, 
with Bh = 0, i.e., 5fi = 0, we find 



ft < Ac , 



a c {T = 0,8n = 0,n)~ i 2*-y/2\eB ± |C(i/2)A* ' r (4.4) 
[0 , Ai>Ac 



where 



- XR\eB ± \ 

Mc — ; ■ (4.5) 

Air - 2 V2RJC(1/2)A« 

These results can be compared with the one in the fully polarized case, 5[M = fi ^ 0. Taking |eBj_| 3> /i 2 and \i 3> a c , 
we find 
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a c (T = 0,5(1 = fx) 



\n\eB A 



4tt~2^BZ\C(1/2)X r ' 
which is half of the value obtained for a c when fj, = Sfi = 0, Eq. ()4.3j) . 



(4.6) 



B. Hall conductivity 

It is not only convenient but also useful to look at possible physical quantities that can be measured in the laboratory 
and give confirmation of the existence of the reentrant and multiple phase transitions we have found above. In most 
realistic experiments, the measured quantities of interest are related, for example, to susceptibilities, like the magnetic 
susceptibility and also quantities related to electric transport. In the presence of a tilted magnetic field, we can have 
each field contributing differently to these quantities. For example, it is expected that the Zceman splitting caused 
by the B\\ component of the field to have important effects in the experiments with graphene and also with other 
planar systems. Based on experimental studies of an in-plane magnetic field in graphene, the Zeeman splitting has 
been shown to be important in both the spin transport and conductance fluctuations properties [l9j . It has also been 
shown that the Zeeman splitting leads to the spectrum of the effective single-particle Hamiltonian exactly as required 
by the observed pattern of quantization of Hall conductivity (2p| . Since the Hall conductivity is a typical quantity of 
importance for analysis when studying these systems, we will here concentrate on the possible effects the reentrant 
and multiple phase transitions we have found can have on it. 

The Hall conductivity, er#, can be defined as [l(| 



where n is the number density, 



\eB ± \ 



n = -^-V e s iR {a c ,T,n,Bj_,5ii) 



(4.7) 



(4.8) 



From the zero temperature limit for the effective potential, Eq. (|4.ip . we find 



a H (o- c ,T = 0,fi,B_i,5fi) 



Ne 2 

4-7T 

Ne 2 

47T 



1 + 2 Int 



sign(Mj.) 



i4- a c 

2\eB±\ 



0(mt - 



1 + 2 Int 



2\eB i I 



ImJ - Wc\) , 



(4.9) 



where sign(x) is the sign function. 
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FIG. 9: The Hall conductivity, at zero temperature, as a function of chemical potential (plots on the left). The vertical thin 
dotted lines separate the regions for which a c changes value (discontinuously) . The plots on the right show the corresponding 
variation of a c with the chemical potential. From top to bottom we have the cases of S/j, = (i.e., B\\ =0), S/j, — 0.7<ro and 
Sfi = 0.9cto. The perpendicular component of the magnetic field is held fixed at the value |ei?x| = Co- 



in Fig. [5] we show some representative cases for the Hall conductivity at T = 0, as a function of the chemical 
potential and for different values for the asymmetry 5fj,. We have again kept fixed, for convenience, the perpendicular 
component of the magnetic field in the value |e£>jj = <Tq as a representative example. The vertical dotted lines in the 
plots on the left indicate the position of the phase transitions, with the shaded areas indicating the chiral symmetry 
breaking regions and the corresponding values assumed for a c - The pattern of chiral phase transition in each case 
is shown in the plots on the right. Note that each time there is a transition, the Hall conductivity can jump by a 
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different factor than in the absence of chiral symmetry breaking. As the chemical potential is increased from zero, 
the largest differences appear at the first transition and at the last transition, when the system enters in the chiral 
symmetric region. This behavior is most different from that seen when computing the Hall conductivity in quantum 
electrodynamics (see, e.g., Ref. pll ]) where there are no gaps and transitions, i.e., where o c = 0. The presence of a gap 
a c and because it can change as the chemical potential is changed, as seen in the previous subsection, is responsible 
for the differences in the Hall conductivity. This is determined by the relations between t^j, and <r c , for a given value 
of B±, as can be inferred from Eq. (|4.9[) . Recalling that fj,^ = fi + Sy and /ij, = fi — Sfx, we find for the Hall conductivity 
the following behaviors as we vary the values of y, and <5/i: 

• When ^ > 0: 



47T 



cr H {(Jc,T = 0,fj,, Bj_,Sfj,) = 



1 + 2 Int 



1 + 2 Int 



^ 2\eB ± \ f 
2\eB, I 



^(Mt-I^c|) 

*(hJ-N) 



(4.10) 



If \jjL\\ > |er c |, we always have that /z-j- > |tr c | and 
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<7H(<r c ,T = 0,h,B±,6/j,) = 2 
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2\eB ± \ 



Int 



%eBA 



(4.11) 



If m > \a c \ and < \a c \, then 



4tt _ f fil — al 

— a H (a c ,T = 0,n,B ± ,6n) = 1 + Int — — 



(4.12) 



— If /if < |cr c |, we always have that \fx^\ < \<r c \, then 



47T 

TVe 2 " 



<jh{<j c , T = 0, /i, Bj_, 5/j,) = 0. 



(4.13) 



When m <0: 
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— If > |(j c |, we always have that /i-j- > |<r c |, then: 



47T 
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aH(<rc,T = 0,/x, Bj_,5fi) = 2 



Int 
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2|eB, I 



(4.15) 



If /i-f- > |(r c | and |/i jj < |(7 C |, we have that: 



47r _ / A 4 ? - CT c 

(Th{(Tc,T = 0,/j,,B±,Sfj.) = 1 + Int I 2 i g I 



TVe 2 



(4.16) 



17 



— If /if < |<7 C |, we always have that \fi±\ < \a c \ and 



—a H (a c ,T = 0,n,B ± ,5n) = 0. (4.17) 

The Eqs. (|4.16p and (|4.17p are identical to Eqs. (|4.12p and (|4.13p . respectively (the sign of /i^ does not interfere in 
these cases). Thus, we have four distinct possibilities for the Hall conductivity, given, respectively, by Eqs. (|4.1ip . 
gH , 6I5])) and (HT71) . 

We also note that by accounting for the effect of thermal fluctuations, we will have two effects on the Hall con- 
ductivity. First, as seen in the previous subsection, it will tend to smooth out the intermediate transitions and only 
the strongest first-order phase transitions tend to remain at a given temperature. Second, including the effect of 
temperature, the Hall conductivity itself is smoothed out. This is shown in Fig. 1101 for the example of asymmetry 
S/j, = 0.9ao that was considered in Fig. [§] Note that for the values of temperature considered, the intermediate 
transition disappears and remains only the first and last transitions, the one entering in the chiral symmetry broken 
region and the one leaving it, respectively. 
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FIG. 10: The Hall conductivity for \eB±\ = erg, S^i = 0.9cr and for T = (dotted line), T = 0.05ct (dashed line) and T = 0.1a 
(solid line). 

V. CONCLUSIONS 

In this work we have studied how an external tilted magnetic field, with nonvanishing parallel and perpendicular 
components to the system's plane, can affect the phase diagram for planar fermions in a GN type of model with 
discrete chiral symmetry. Though in this work we have concentrated on the transition patterns found for the chiral 
fermionic condensate, we can easily extend our work to include also a superconducting gap (22j . We will investigated 
this in details elsewhere. 

As far the chiral symmetry is concerned, the parallel component of the magnetic field leads to a Zeeman effect in 
these planar systems that tends to weaken the chiral symmetry. Thus, the chiral symmetry phase transition, which in 
condensed matter systems can be associated with an insulating to metal type of transition, may happen at a smaller 
critical temperature in the presence of a Zeeman field. This behavior, due to an increasing magnetic Zeeman field, is 
exactly the opposite to what is observed when a perpendicular magnetic field is applied to the system (see e.g. [6|,|7|). 
In a perpendicular magnetic field the chiral symmetry breaking becomes stronger as a consequence of the effect of the 
magnetic field, a behavior associated with the magnetic catalysis in the system (note that magnetic catalysis is not 
an effect exclusive for fermionic systems, but can also occur in bosonic systems as well (23|, enlarging the symmetry 
broken region in those systems). Thus, in the presence of a perpendicular magnetic field the transition happens at a 
higher critical temperature. We have shown that this competing effect of each component of the external magnetic 
field applied to the system can produce a rich phase diagram, with the production of multiple critical points and 
reentrant phase transitions. 
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Note that there are many systems that can display multiple phase transitions, most particularly in low energy 
condensed matter systems, where they can be found more easily. At the quantum field theory level, reentrant phases 
have been shown to be possible in systems with coupled scalar fields in the nonrelativistic limit [24| . Multiple phase 
transitions can in principle also be found in relativistic systems. This can be the case, for example, in systems with 
two or more coupled scalar fields (see, for example, [25[ and references therein). Multiple critical points in effective 
models for quantum chromo dynamics have also been shown to be possible in Ref. (26|. However, to our knowledge, no 
reentrant phase transitions of the type we have found here have been shown before to appear in fcrmionic relativistic 
systems. 

As a way to determine the presence of multiple critical points and reentrant phase transitions, we have computed the 
Hall conductivity. We have seen that whenever a first-order phase transition happens, the Hall conductivity jumps by 
a different amount at the critical point. This is a consequence of the different value assumed (discontinuously) for the 
chiral vacuum expectation value tr c , as analyzed in subsection IIV Bl Since the Hall conductivity an is an important 
measurable quantity for diverse condensed matter systems, like for semiconductor materials, superconductor films 
and graphene, our results indicate a way of how the presence of a gap in the system and quantum phase transitions 
may reflect on the measurements of an- It is tempting to try to directly apply our results to the physics of graphene, 
for instance. In these systems we can have both components of the magnetic field acting concomitantly. However, 
in most studies (sec, e.g., [13) the Zeeman energy term (which for us here is given by the asymmetry term S/i) 
is much smaller than the energy splitting between Landau levels, ^/2|ei3j_ | . This is easily understood if we use 
typical values found, for example, for graphene. In this case, 5/i = gfiBB\\/2, with g = 2 and hb = e/(2m e ), where 
m e = 0.511 MeV is the electron rest mass. Using also that e ~ 1/V137 and the conversion factor 1 tesla ~ 692.4 eV 2 , 
we find 5/j, ~ 0.058 Bn [T] meV, where B[T] is the absolute value of the magnetic field in teslas. For the Landau energy 

level splitting between the zeroth and first levels, we find that Ell = \j2v F \eB^\ ~ 36.3^/-Bj_ [T] meV, where we have 
re-established the Fermi velocity in the expression and used that vf ~ c/300 for graphene. We then have that for 
typical laboratory magnetic fields of (9(10) teslas and for a not too much tilted magnetic field (such that B± ~ B\\), we 
find that 5fj, <C Ell- So, indeed, under these circumstances, the Zeeman energy term becomes negligible. However, 
in our case, we can find already reentrant and multiple phase transitions for asymmetries as small as Sfi ~ O.loo and 
at a temperature T < O.Olero- Taking a perpendicular component for the magnetic field of vf^e-BjJ ~ 10c 2 (where we 
have re-established again the Fermi velocity in the expression, as appropriate for fcrmionic condensed matter systems) 
and for a typical gap energy <jq of around 10 meV [28[, we find T < 0.1 K and a required field with By ~ 17.3 teslas 
and B±_ ~ 1.5 teslas. This corresponds to a highly tilted magnetic field in the direction of the system's plane (see 
Fig. [1]), with tan0 ~ 8.7 x 10~ 2 . However, these estimates for the required values of magnetic fields are reasonable 
and low enough to be achieved in the laboratory and this does not seem such a highly special situation to be produced 
in practice. Besides, with the increasing precision being reached in the most recent experiments, it may be possible to 
realize such an experiment with superconductors films under these tuned conditions and probe the type of transitions 
we have found in this paper, for example through the measurements of the Hall conductivity, as we have suggested. 

In conclusion, we have seen that the opposite effects caused by magnetic fields, when applied parallel and perpen- 
dicular to the system's plane, can lead to a rich pattern of phase transitions. The results we have found here can have 
immediate applications in the context of the physics of planar condensed matter systems, including graphene, planar 
films of conducting polyacetylene and high-temperature superconductor films. In all these systems, the viability of 
their application as modern semiconductor devices is directly related to the ability of producing and controlling a 
gap. Our results indicate that this control can be achieved through an appropriate doping of the system and under 
parallel and perpendicular magnetic fields, that are applied simultaneously and independently, when properly tuned. 
If these conditions can be achieved, then our results can open interesting possibilities for the uses of these type of 
materials in practical applications as novel electronic devices. 
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